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^ ' Abstract 

Q^' The parton model as developed for semileptonic B decays is applied to the in- 

clusive decay B J/tpX. We calculate the momentum spectrum of the J/^ using 
a one-parameter distribution function for the heavy quark, taken from production 
^ ■ experiments, and compare our results with recent data from CLEO which fixes the 

' distribution parameter Sp. An analogous calculation is carried out in the ACCMM 

model where the data determines the Fermi motion parameter pf. 

The models give a good description of the data, provided in each of them the 
parameter in the distribution function is chosen to correspond to a softer b quark 
momentum distribution than that commonly used in studies of semileptonic B decays. 
In particular we arrive at Sp = 0(0.008) and pj = 0(0.55 GeV), respectively. The 
latter value {pf = 0.5 GeV) produces in the ACCMM model 10^ x \Vub/Vcb\'^ = 1.03, 
which removes a former supposed discrepancy with the ISGW model. Finally, the 
strength of the effective color singlet coefficient is found to be \a2\ = 0(0.28) in both 
models. 
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1 Introduction 



The distribution of charmonium states in B meson decays is of special interest, because it 
provides a testing ground for color suppression in B decays already attainable in the present 
experiments. At the b mass scale, a leading order calculation generates a consistent pattern 
of color suppression for the production of S-wave charmonium states. In the context of 
the 'color singlet mechanism' the production is related to the decay of a 6 quark in the B 
meson, at short distances, into a color singlet cc pair plus other quarks and gluons. The c 
and c have almost equal momenta and reside in the appropriate angular-momentum state. 
This subprocess is dominated by the color suppressed internal W -exchange diagram, where 
hard gluon exchanges lead to an effective neutral current Deviations from this leading 
order 'singlet mechanism', which arise from relativistic corrections and soft gluon induced 
fragmentation of the cc into charmonium, affect the normalization but not the structure 
of the amplitude b ^ J/ip X {in which X sums over hadronic states) [0. 

Several authors made predictions for the branching ratio of the direct decay b —* J/ip X 
to leading order in as [1-5]. Recently a next-to-leading order QCD calculation was reported 
p. The authors pointed out that the process under consideration cannot be explained 
presently by a standard application of perturbative QCD. The difficulties entering the 
analysis arise from the strong suppression of the leading order color singlet Wilson coeffi- 
cient, causing considerable cancellations in three different orders of as- Alternatively one 
may use a phenomenological factorization approach and take into account deviations from 
this prescription by replacing the color singlet renormalization coefficient by an effective 
neutral current coefficient 02, which has to be determined from experimental data in a 
model dependent way. 

Both ARGUS and CLEO reported on inclusive B decays to J/ip, where they identified a 
sizable component of decays with three or more particles in the final state J^, . In a recent 
publication CLEO presented an analysis based on a data sample, which is one order of 
magnitude larger than those of previous studies, corresponding to an error reduction by a 
factor of 2.4 . A comparison between theory and experiment requires the branching ratio 
resulting from direct decays. Therefore the group corrected the measured J/ip spectrum 
for the 'feed-down' modes B iIj{2s) X and B — > Xci^. As a result they found the 
direct branching ratio (0.80 ± 0.08)%. A theoretical analysis carried out by the CLEO 
group following Ref. yields 0.75%, in good agreement with the measurement, where the 
phenomenological factorization prescription has been applied using a value |a2| = 0.23 . 

The latter value has to be modified when incorporating bound state effects for the initial 
B meson, which have to be included in order to gain a theoretical prediction of the J/ ip 
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momentum spectrum within an inclusive approach. So far there is not available a detailed 
fit to experimental data of momentum spectra obtained from theory. Palmer and Stech 
have made a preliminary attempt using a simple wave function formalism |T^. In this 
paper we will investigate two different approaches in order to account for the experimental 
results. 

As stated above a sizable component of the momentum spectrum is due to nonresonant 
multi-particle final states. Consequently, for a wide range of phase-space an inclusive 
description using quark-hadron duality, which was extensively applied in the inclusive 
semileptonic decays of B mesons, will be appropriate. 

The semileptonic decays involve large energies and momentum transfers of the weak 
current over most of the phase-space. Several groups noticed that these kinematic regions 
investigate the light-cone behaviour of currents for which the methods of deep inelastic 
scattering are valid. There are various approaches available which derived useful results 
describing the lepton spectrum in the decay B Xu(c)ez/. These approaches use different 
formalisms, but they still follow two basic steps. 

(i) The decay of the b u{c) which are propagating as free particles is considered and 
the resulting spectrum is folded with the momentum distribution of the b quark obtaining 



moments of the quark distribution 12]. A similar approach is given by the ACCMM 



model [T^, in which one accounts for the binding effects by treating the spectator quark 



as particle on the mass shell with nonzero momentum and averages over this momentum. 

(ii) After the introduction of the point-like behaviour one still has to calculate the expec- 
tation value of the bilocal transition operator between the initial hadronic state to gain 
the width of the decay. In this context the new development is the application of the 
operator product expansion (OPE), which involves an expansion in inverse powers of the 
heavy quark mass rrih incorporating the formalism of the Heavy Quark Effective Theory 
(HQET) Ijl^. The calculation of the lepton energy spectra requires a modified expansion 
in powers of 1/(1 — y)mh with y being the normalized lepton energy [18-23]. 

The analytic form of the structure function of the heavy quark within the B meson controls 
the endpoint behaviour of the leptonic spectrum, since in this region the OPE does not 
converge. It has been shown that the prediction of HQET far from the endpoint gives 
approximately the same shape of the spectrum as the ACCMM model with a Fermi mo- 
tion parameter pf ~ 0.3 GeV [^]. Yet the model dependence of this region is quite small 



25| . Therefore it is of great relevance to test the approaches with decay channels of the B 



meson different from the semileptonic one in order to extract the momentum distribution 
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function of the heavy quark in a direct way and then compare it with the distribution 
function obtained in semileptonic decays. 

The inclusive decay B J/ipX provides a well suited testing ground since, as we 
shall show, in this process the dependence on the structure function of the h quark is more 
direct, i.e., the J/ip momentum spectrum is proportional to the distribution function. 

A QCD based analysis of the decay using HQET is only of limited validity, since due to 
the smaller energy release the convergence of the OPE is less reliable than for semileptonic 
decays. Therefore one has to resort to a formalism in which the structure function can be 
modeled in terms of parameters that may be obtained from experiment. We will apply 
two different approaches. In Section ^ we present a field-theoretical version of the inclusive 
parton model, which two of us together with Jin already applied for semileptonic decays 



n, |T^. In Section ^ we calculate the J/^ momentum spectrum within the framework of 
the ACCMM model. In each of the analyses we determine the distribution parameter of the 
structure function from a comparison to the recent CLEO data. The overall normalization 
of the theoretical spectra further provides information respecting the value of the effective 
color singlet coefficient \a2\. Finally we use our results to extract the value of |V^fe/V^b| in 
the ACCMM model from the semileptonic decay channel of the B. This procedure requires 
some additional remarks which are given in Section ^. In Section ^ we compare our results 
with experimental data for the polarization of the J/ip. The summary can be found in 
Section 



2 B ^ J 1^ X in the Inclusive Parton Model (IPM) 

2.1 Calculation of the Differential Branching Ratio 

In order to investigate the interaction which induces charmonium production in B decays, 
we start at the quark level at a high-energy scale where the Hamiltonian is known. This 
is renormalized to lower energies to produce an effective four-fermion interaction 
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{s-^d}, (1) 



7^ = 7^(1 — 75), where operators arising from penguin and box diagrams have been ne- 
glected. The first term is a color singlet, the second a color octet operator. The renormal- 
ization functions (Wilson coefficients) Cj(/i) have been computed up to the next-to- leading 
order corrections in Ref. p6| . 

The effective interaction (|lD is evaluated using a factorization prescription by which the 
amplitude can be written as a product of matrix elements of current operators. Deviations 
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from this prescription are parametrized by substituting the color suppressed singlet coef- 
ficient by a free parameter 02 which is written in terms of the effective number of colors 
(1/.^) and has to be determined by mo del- dependent fits to the experimental data, 

C2 + j a2 = C2+ ^Ci . (2) 

02 is equivalent to the coefficient introduced for type-II processes in the factorization model 
of BSW 1^ for exclusive decays. Using CLEO data from two-body decay modes with J/ ip 
mesons in the final state, an analysis within the BSW model yields 



lasi = 0.26 ±0.01 ±0.01 ±0.02, (3) 

where the second systematic error is due to the B meson production fractions and life- 
times. Theoretical uncertainties mainly due to hadronic form factors are not included. A 
corresponding analysis using the model of Neuhert et al. provides a central value 
\a2\ = 0.23, whereas that of Deandrea et al. yields |a2| = 0.25 (for \Vcb\ = 0.041, 
tb = 1.44 ps and foiD*) = 220 MeV). 

In this paper we extend the factorization model for exclusive two-body decays to an 
inclusive picture at the quark level. This is reasonable because, as was stated in Ref. p2 



many-body final states will most likely start as two color singlet quark antiquark pairs, in- 
cluding intermediate massive resonances, where strong phases from final state interactions 



disappear in the sum of all states. The authors of Ref. |32| made use of this inclusive pic- 
ture to determine the lifetime ratio t{B^)/t{B^). As we intend to reproduce not only the 
branching ratio but also the momentum spectrum for the decay B ^ J/ipX, we have to 
incorporate the momentum distribution of the b quark in the B meson. In our analysis we 
will let the effective color singlet coefficient \a2\ to fioat within the range which is covered 
by the various models for exclusive modes. 

We consider the matrix element of the Hamiltonian (|I]) for our process (in the explicit 
notation we only refer to the dominant decay involving an s quark in the final state). 



(J/^X,|K//|i?) = a2^K6K:(W|c7M0)(X,|s7^6|S) 



(4) 



where we make use of the phenomenological factorization prescription. The second term 
of Eq. (|1|) does not contribute between color singlet states. 

For the matrix element of the cc{lS) at the origin we use the current- identity for pure 
vector-like states, 

(0|c(0)7M0)|JM=£^/v,M^. (5) 
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This identification is valid, because the time-scale of the interaction to be considered is 
larger than the scale for the formation of the singlet state {tint > '^/M^). is the J/ip decay 
constant and can, to leading order in the relative heavy quark velocity v, be related to the 
radial part of its nonrelativistic wave function. In principle one can determine from the 
electromagnetic decay J/ip ^ e^e~ , but practically one is confronted with the ignorance 
of higher order QCD corrections to the decay rate which might be of great relevance in 
view of existing large leading order corrections. Within the studies of exclusive two-body 
decays B — * J/ipM the central value of the decay constant was fixed at = 0.384 
GeV which results from the omission of any QCD correction (as of relativistic corrections 
to the wave function of the J/ip)- This analysis refers to the 1990 Particle Data Group 
value for the branching ratio of J/ip ^ e~^e~ . Therefore when extracting \a2\ from inclusive 
B decays we will use the above value in order to compare the result with Eq. (|^). This 
procedure requires some additional remarks which concern both the studies of exclusive 
and inclusive decays. 

Relating the parameter which contains the nonperturbative effects in the production of 
IS'-charmonium with the J/ip decay constant implies the neglect of color octet contribu- 
tions to the fragmentation process in the explicit calculation. In Ref. it was argued in the 
framework of nonrelativistic QCD (NRQCD), that the fragmentation of a {cc)^{^Si) state 
into a J/ip in the long-distance scale may give sizable contributions to the B J/ipX 
decay rate. This feature occurs, despite the fact that the color octet operator is of the 
order (9(t>^) respecting the velocity counting rules in the NRQCD, because of the large 
suppression of the color singlet Wilson coefficient. 

As the structure of the amplitude for the effective point-like decay b ^ J/ipX re- 
mains unchanged, considering nonleading effects in the fragmentation process only affects 
the normalization of the B J/ipX decay spectrum. Using the leading order value 

= 0.384 GeV therefore means that the effective neutral current coefficient \a2\, which 
is determined from a comparison with experimental data, includes contributions due to 
the charmonium fragmentation. Consequently, since these contributions are different for 
various charmonium states [^, the universality of |a2| is limited to B decays involving fixed 
S-wave resonances. 

Reducing the matrix element (^) of the mesonic decay B ^ J /ip Xg to that of the point- 
like effective neutral current subprocess h J/ips (see Fig. Ic) in the valence quark 
approximation one obtains 

{J/i,X,\neff\B) = a2V,,V:,M^U^e;sY{l - l,)h . (6) 
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This expression for a free quark decay has to be modified by incorporating bound state 
corrections due to the strong interaction between heavy and hght quark in the initial meson. 
In the intuitive parton model one has to take the incoherent sum over all subprocesses in 
which a distribution function f{x) is introduced as a weight, 

Tb = dxf{x)Tb{x), where / f{x)dx = l. (7) 
Jo Jo 

In Eq. d^) transverse momenta of the heavy quark relative to the B meson are neglected 

which leads to the Lorentz invariant prescription p^^ = xPq. This approach corresponds to 

an equal velocity approximation, i.e., valence quarks and B meson have the same velocity, 

and is valid at the infinite momentum frame. We therefore proceed as usual computing a 

Lorentz invariant quantity and then use it in any other frame. 

The basis for our approach is given by a field-theoretical version of the parton model. 

If we square the second matrix element on the right-hand side of Eq. (^) and sum over all 

final states, which guarantees incoherence, we produce the hadronic tensor corresponding 

to the effective transition B — > Xj, 



W,. = -^J dSe"^y{B\ [j,{y), jt(0)] \B) , 



where jfj.{x) =:qf'^^h{x) : is the left-handed effective neutral current. The same tensor 
structure we encountered in the semileptonic B meson decays [11-14]. Substituting for the 
leptonic tensor the corresponding expression of the J/ip current, 

= 2AV.f\^\a,\^flMl ^-g,^ + , (9) 

we can write the differential rate for the decay B Xj J /ip in the restframe of the B in 
exact analogy to the semileptonic decay, 

<'r,«-.,.«^f|^V»-5. (10) 

The general structure of the hadronic tensor reads 

= _gf^^W, + ^P^BPBW2-^e^^"P^PBo.qpW, (11) 

B B 
B B B 

where in our case q = k^. Introducing the light-cone dominance as in Ref. [11-15] allows 
to refer the tensor corresponding to the transition B —>■ Xf in the decay B Xf J /ip to 
the distribution function f{x) of the heavy quark momentum, 

W^,j = A{S^p,x-ie^p,x) J dxf{x)P^{xPB-k^Ye[{xPB-k^)omxPB-k^f-m)], (12) 
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with 



fipi/X 



Qppgux - g^lugpx + g^xQup and e{x) 



+ 1. 
-1 



X > 
X < . 



(13) 



From Eqs. ([TT| ) and ([T2|) one obtains the structure functions Wi of the hadronic tensor in 
the restframe of the B meson, which two of us together with Jin derived in Ref. [1TI| , |l^ in 
exact analogy for the semileptonic decay rate, 
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W, = 2[/(x+) + /(x_)] 
-4 



x+ — x_ 
where we have defined 



[/(x+) - /(x_)] 



W2 



H^6 = 0, 

(|Ps|=0) / 



(14) 



E^±\\k^\^ + m) . (15) 



For the case of a massless final state quark x± are identical to the usual light-cone variables. 
The dependence of the distribution function on the single scaling variable x is a consequence 
of the light-cone dominance, since in this framework the structure function /(x) is obtained 
as the Fourier transform of the reduced bilocal matrix element, which contains the long- 
distance contributions to the hadronic tensor 1 11 311, 



/(x) = / d{yPB)e 



ix{yPB 



(16) 



The terms proportional to /(x_) in Eq. (14) are a result of the field-theoretical approach. 
The kinematical range for x_ belongs to a final state quark with negative energy. Therefore 
the corresponding terms can be associated, formally, with quark pair-creation in the B 
meson (see Fig. lb), whereas the dominant terms proportional to f{x^) reflect the direct 
decay of Fig. la. 

Including the small f{xJ) term as well as the CKM-suppressed transition c — > (i the 
Lorentz invariant width of the decay B ^ J/ip X reads 

|2 



Eb ■ dT(B. 



(B^.J/i,X) 



E 

f=s,c 



\c 



271^ 



dx /(x) 



Pb{xPb - 



— {Psk^){xPBk^-M^) 



xe[{xPB - k^)o]S^^^ UxPb - k^f - m) 



d^k^ 
2E,i, 



(17) 



with 



\Cf 



G p 
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\Vcb\\Vcf\M.^f^\a2\ 



8 



Here the modification from tlie field-tlieoretical approacli simply enters in form of the step- 
function e{x) which exactly provides the additional terms f{xS) in the tensor coefficients 
Evaluating Eq. (P!7| ) in the restframe of the B meson we arrive at the formula for the 
J/ip momentum spectrum, 



31^1 + W2 



2 



Ml 



(19) 



within the framework of the inclusive parton model. Using Eq. the corresponding 
differential branching ratio can be written as 

1 '^^^--^l^-)(\k^l\f;,\ = Q)= (20) 



d\k 



^1 



2 17.12 



f-, 271 Mb 
within the kinematical range 



< < ^^^^^ , b^i^-rrij.. [ZL) 



B 



To compare our result ( pOf ) with data from CLEO a Lorentz boost has in addition to be 
performed from the restframe of the B meson to a -B produced at the T(4S') resonance 
{\pb\ = 0.34 GeV). The explicit form of the boost integral we give in Section in the 
context of the ACCMM model (see Eq. (H)). 



All terms are known in this decay except for the product {\a2\f1p) and the structure function 
which occurs with two arguments. We note that when neglecting the small f{xJ) contribu- 
tion the decay spectrum is directly proportional to /(a;+). Thus measuring the momentum 
spectrum of the J / we can read from the data the distribution function. The extraction 
of /(x+) is direct, in contrast to the extraction from the electron energy spectrum in the 
semileptonic decay of the B meson which involves an integral over the structure function. 

In the following Section we compare the predictions of Eq. (|20| ) with the existing ex- 
perimental data. 



2.2 Analysis and Numerical Evaluation 

The probability distribution function (0) is not Lorentz invariant and is defined in the 
infinite momentum frame. It cannot be transformed to the restframe of the B meson (or 
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a frame which corresponds to a 5 meson produced at the T(4s) resonance), because it 
involves an infinite sum of quark-antiquark pairs whose calculation requires a complete 
solution of the field-theory. In the absence of direct measurements of the distribution func- 
tion we use a one-parameter Ansatz and fix the distribution parameter by comparing our 
results with data. Referring to theoretical studies which pointed out that the distribution 
and fragmentation function of heavy quarks peak at large values of x ||3^, ^ , as a working 
hypothesis we assume that the functional form of both is similar. The latter is known from 
experiment and we shall use the Peterson functional form [35-37] 

with Ep being the free parameter and A^^ the corresponding normalization constant. This 



form we already applied in the semileptonic decays of the B meson ||TT] , |T2|. 

In Fig. 2a we show the distribution function for various values of Ep. The kinemat- 
ical range for the two arguments x±{mf) of the distribution function according to Eq. ( ^11) 
reads 

— ^— — -<x+<l, — ^— ^ — ^ < x_ < — ^— — -. (23) 
Mb ~ Ml - - Mb 

As pointed out in Fig. 3 the variable x_ only occurs with values at which /e(x_) is small. 

Therefore the corresponding contribution to the differential branching ratio ( pOD is small. 

We use the distribution ( p2l) to fit the measured momentum spectrum of the J /if) which 
was given by the CLEO group 0. As argued above we fix = 0.384 GeV and vary \a2\ 
according to the range covered by the various models for exclusive decays. 

The shape of the theoretical spectrum is determined by the distribution parameter Ep, 
which is illustrated in Fig. 4 where we show the spectrum for \a2\ = 0.275 and several values 
of Ep. A general feature of the analysis is our difficulty in reproducing the data over the 
whole range of phase-space. Confronted with this problem, we lay greater emphasis on the 
appropriate description of the /ow momentum range {\k^\ < 1.4 GeV) in our simultaneous 
fits of \a2\, which in this decay has the meaning of a pure normalization constant, and 
the distribution parameter Ep. Within this region the J/ip spectrum obtains a sizable 
contribution from decay channels containing three or more particles in the final state 
(where higher K* resonances are assumed to be unimportant), whereas the high momentum 
range is purely determined by the exclusive two-body decays B ^ J/if) K^*^ . Therefore 
incoherence as a necessary ingredient of the IPM is limited to the former region. 

Furthermore we demand the reproduction of the measured branching ratio for the 
decay under consideration. Thus if we apply small values for the distribution parameter 
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Sp < 0.006 we cannot account for the low momentum region, which is underestimated 
significantly within the theoretical spectrum (see Fig. 5). On the other hand, demanding 
the reproduction of the branching ratio, a distribution corresponding to Sp > 0.010 would 
imply a value \a2\ > 0.30, which is beyond the range given from the study of exclusive 
two-body decays. Therefore, in spite of the difficulty to reproduce the data accurately 
over the whole range of phase-space, the investigation of the decay B ^ J/ipX is more 
restrictive with regard to the distribution parameter than the semileptonic decay of the B 
meson which involves an integral over the structure function. 

Considering the fact that the IPM does not include hadronization effects and conse- 
quently yields an averaged spectrum, a satisfactory description of the dataset is achieved 
for Sp = 0(0.008), \a2\ = 0(0.285) (compare Fig. 5) when using current masses for the 
final state quarks. Applying constituent masses would involve an incorrect position of the 
maximum of the momentum spectrum and, in addition, a sizable underestimate of the 
high momentum range (see Fig. 6). Yet also for Sp = 0.008 there is apparent a moderate 
systematic underestimate of the latter region. Since the model presumes the existence of 
multiple final states it may not hold for the high momentum range {\k^\ > 1.4 GcV) where 
two-body decays determine the spectrum. Following the CLEO analysis, as a trial we 
subtract the expected contribution of the exclusive decays B K^*\j/%p from the data for 
the semi- inclusive momentum spectrum. Carrying out the fit for this modified momentum 
distribution we observe that the agreement between theory and experiment is improved for 
large values of the distribution parameter Sp > 0.012 (however only when applying con- 
stituent masses for the final state quarks) , corresponding to a soft momentum distribution 
of the heavy quark. In Fig. 7 we show the modified spectrum for |a2| = 0.275 and £p as 
stated there. This additional fit may serve as an indication that a value Sp = C(0.014) 
has to be suggested within the model when restricting the analysis of decay spectra to the 
'incoherence region' where multi-particle final states exist. Then additional information 
from a model describing the dominant two-body decay channels has to be implemented to 
reproduce the data over the whole range of phase-space. 

Finally, the errors in the data, although substantially improved, are still significant and 
a crucial test will be possible, when the error bars are further reduced. It is of special 

— * 

interest to establish whether the deviations at \k^\ ~ 0.5 GeV from the smooth shape 
of the spectrum, which is predicted in our inclusive approach, survive and to extend the 

— * 

analysis of the composition, in terms of resonances and continuum, for \k^\ > 1.4 GeV. 
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3 B ^ J/^Ij X in the ACCMM Model 

3.1 Calculation of the Differential Branching Ratio 

A second approach which allows us to determine the momentum distribution of the J/ip 
in the semi-inclusive decay of the B meson is given through the ACCMM model [rB]. In 
this model the bound state corrections to the simple quark picture are incorporated by 
attributing to the spectator quark a Fermi motion within the meson. The momentum 
spectrum of the J/ ip is then obtained by folding the Fermi motion with the spectrum from 



the b quark decay. In Ref. the shape of the J/ijj momentum distribution resulting 
from Fermi momentum smearing has been given in the restframe of the B meson, but 
without a comparison with data, and without consideration of the light spectator quark 
mass. Moreover the authors applied the constituent mass for the strange quark in the 
final state, whereas we will show that a satisfactory reproduction of the data can only be 
obtained applying the current mass. 

According to Eq. the J/ ip is treated as a color singlet current using the factorization 
assumption of Eq. (^). The spectator quark is handled as on-shell particle with definite 
mass rrisp and momentum \p\ = p. Consequently, the b quark is considered to be off- 
shell with a virtual mass W given in the restframe of the B meson by energy-momentum 
conservation as 

W'ip) = M| + m% - 2MB^m% + p^. (24) 

Altarelli et al. introduced in the model a gaussian probability distribution (pip) for the 
spectator (and thus for the heavy quark) momentum, 

0(p) = -^exp(V/p;) , (25) 

normalized according to 

/ dpp'^(j){p) = 1 . (26) 
Jo 

Here a free parameter p/ is adopted for the gaussian width which has to be determined by 
experiment . 

The main difference between the inclusive parton model and ACCMM is that the latter 
one must consider a b quark in flight. We therefore start from the momentum spectrum 
of the J/tp resulting from the decay b —>■ qf J/ip [f = s, d) of a b quark of mass W and 
momentum p which is given by 
jp(/) 



d\k^\ ^ ^ ^ kf'>{p)-\k'^''\p)\ 



(27) 
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Here we have defined 



Fq is the width of the analogous decay in the restframe of the heavy quark where we have 
confined ourselves to the leading order result, 



9{x) 



1, 
0, 



X > 
X < . 



f28) 



2 u(f) 



\C£kl_ 
27r 



m 



+ (w'-mj-Ml) 2 + 



(29) 



with \Cf\ according to Eq. (|18|) and with the momentum k^^^ of the J/i/j 

1 

1 



if) 



2W 



if) 



(30) 



For a vanishing mass mj of the final quark Eq. (|29|) is equivalent to Eq. (5) in Ref. 0. 

In Eq. ( p7|) A;±' give the limits of the momentum range which results from the Lorentz 
boost from the restframe of the h quark to a frame where the h has a nonvanishing mo- 
mentum p, 

kl'f\p) = ^{E,k\l^±p4f^) , (31) 
and is the corresponding Lorentz factor, 



7b 



W 



Eh 



(32) 



To calculate the momentum spectrum of the J/ip from the semi-inclusive decay of the 
B meson one has to fold the heavy quark momentum probability distribution with the 
spectrum (^) resulting from the quark subprocess. Performing this we finally arrive at 
the expression for the differential branching ratio for a B meson in flight. 



1 dT 



d\k^\ 



[\k^\, \pb\ 



(33) 



d\k'^ 



Fmax 



dp p^(j)(p) — IT— 

k4\k'^\) - \k.{\k'^\)\ { 'd\k^\ 



dr 



if) 



Here p^l^ is the maximum kinematically allowed value of the quark momentum p, i.e., 
that which makes W in Eq. (^^ equal to W = rrif + M^, 



pif) 

fmax 



2M 



B 



{Ml + m% - {mf + M.^ff - Am%M^ 



(34) 
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The first integration in Eq. (^) results from tlie transformation from the spectrum for a 
B meson at rest to the spectrum for a B meson in flight, where 

k±{\h\) = ^ {Eslhl ± \Pb\E^) , k^-i\\t^\) = minik^lQ), kill} , (35) 



kmax being the maximum value of the J / if) momentum from the decay B ^ J/'ip Xj in the 
restframe of the B, 



2M. 



B 



{Ml + Ml- 4i) - AMlMl " , = (^/ + m,,f . (36) 



In the following Section we make use of Eq. (|33D to compare the model predictions with 
experimental data. 

3.2 Analysis and Numerical Evaluation 

Both models, the inclusive parton model as well as the ACCMM model incorporate the 
bound state structure of the B meson by postulating a momentum distribution for the 
heavy quark. Introducing in the latter model another x- variable as the ratio x = W/Mb, 
the appropriate distribution function for the relative mass fraction w{x) of the b quark in 
the restframe of the B meson is given as 



p{x) 



2 

where the corresponding normalization reads 



2 4 ■ 

TP Tft 



(37) 



1— r?isp/Ms 

w{x)dx = l. (38) 



In Fig. 2b we show the distribution function w{x) for m^p = 0.2 GeV and various values 
oi Pf. From Eqs. ( |H7| ) one realizes that for given on-shell masses the Fermi parameter pf 
determines the average value (W) = Mb{x) of the virtual b quark mass and hence also the 
total decay width. 

Both models have the advantage of avoiding the mass of the heavy quark as an in- 
dependent parameter. As a consequence, the phase-space is treated correctly by means 
of using mesonic degrees of freedom. Nevertheless the endpoint of the phase-space in the 
ACCMM model is slightly different from that in the IPM because of the on-shell mass 
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nisp of the spectator quark from which the minimal invariant mass square of the hadronic 
system S^^^ results as given in Eq. (^). 

The shape of the momentum spectrum in the decay B J/ipX is mainly determined 
by the value of p/. The extraction of this value from a comparison with experimental 
data is important not only for explaining the decay itself and thus testing the factorization 
assumption of Eq. (||), but also for the determination of iKifc/Kb] from the endpoint region 



of the inclusive semileptonic B — > Xc(m)/z/ decay spectrum. As recently stressed in Ref. p9 



the experimental extraction of from semileptonic decays has been ambiguous, because 
various parameters of the model were fitted simultaneously to the lepton energy spectrum 
where in addition the perturbative QCD corrections are important in the endpoint region 
(especially for h ^ u). Furthermore, as stated in Ref. [25|, for a large range of the phase- 
space the bound state corrections are of minor importance in this decay channel. 



Whereas usually pf = 0.3 GeV is used for the data analysis the authors of Ref. ^9 
calculate the Fermi parameter theoretically in the relativistic quark model and obtain 
Pf = 0.54 GeV. A CLEO analysis |0 of the endpoint lepton spectrum in the semileptonic 
decay channel employing the common value pf = 0.3 GeV (and m^p = 0.15 GeV) yields 
a discrepancy between the ACCMM and the ISGW model of Isgur et al. in the 

determination of 

lO^xlKfe/KbP = 0.57 ±0.11 (ACCMM) 

= 1.02 ± 0.20 {Isgur et al.) . (39) 



Using however pf = 0.5 GeV the authors of Ref. |^ arrive at the value 10 x \Vub/Vch\ = 
1.03 within the ACCMM model, obtaining a good agreement with the ISGW model. 

The decay B —>■ J/ip X provides an independent way of extracting the Fermi parameter 
Pf. In contrast to the lepton spectrum in the semileptonic decay the shape of the J/ip mo- 
mentum spectrum is highly sensitive to this parameter over the whole phase-space. This 
is analogous to the determination of the parameter Sp in Section P]2 . 



Employing the spectator distribution function ( p5|) we calculated the momentum spectrum 
for the decay of B mesons produced at the T(4s) resonance. Fig. 8 shows the comparison 
with the CLEO data for = 0.384 GeV, m,p = 0.2 GeV, m, = 0.125 GeV, = and 
various sets of values for pf and \a2\ as stated there. Here we want to emphasize that 
analogous to the IPM the best fit is achieved when employing current masses for the final 
state quarks. This is a consequence of the fact that within the theory these are treated 
as free particles without consideration of fragmentation. The application of constituent 
masses would yield a sizable underestimate of the high momentum region (see Fig. 9). 
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In addition to the final state quark masses mj the spectator mass also determines the 
range of phase-space and the position of the maximum (compare Fig. 10). The latter is in 
accordance with the data for rrisp = 0.2 GeV. 

In Fig. 8 one can see that the agreement of the theoretical spectrum with the data 
is very good provided we choose the value pf = (9(0.55 GeV). To demonstrate the high 
sensitivity of the spectrum on this parameter over the whole range of phase-space we 
present in Fig. 11 spectra for different values of pf with and \a2\ fixed. It is obvious that 
the shape of the experimental J/ip momentum distribution (compare Fig. 8) cannot be 
reproduced in the model when using values for the Fermi motion parameter significantly 
smaller than pf = 0.5 GeV. Especially the sizable contribution in the low momentum range 
requires a soft probability distribution of the heavy quark momentum. 

As a result, we conclude that the value pf = (9(0.55 GeV) for the Fermi motion parame- 
ter of the ACCMM model is highly favoured when applying this model to the semi-inclusive 
decay B —>■ J/tp X . This value is in exact accordance with the one calculated recently in 
Ref. 1^ from the relativistic quark model. As mentioned above the discrepancy between 
the ACCMM model and the ISGW model concerning the determination of |V^6/V^fe| dis- 
appears when using our favoured value for pf in the study of the electron spectrum from 
inclusive semileptonic B decays. 

A final feature of our investigations is related to the determination of the effective 
color singlet coefficient \a2\. Using = 0.384 GeV, tb = 1-54 ps, and \Vcb\ = 0.043, we 
obtain jasl = (9(0.28) for the best fit of the ACCMM model to the CLEO data. This is in 
agreement to the value of \a2\ which we found in the IPM. The study of an inclusive decay 
compared to exclusive decays in this context provides the advantage of solely involving an 
analysis that is independent of mesonic form-factors. 



4 Remarks on the Extraction of 

The application of the value for pf, which we obtain from our investigation of the decay 
B — > J/ipX, to the semileptonic decay channel requires some additional remarks. As 
pointed out in Ref. |^ the universality of the momentum distribution function describing 
the motion of the b quark in the B meson only holds for different decay processes when 
referring to a fixed final state quark mass. 

The common procedure, followed in the analysis of the semileptonic decay B Xlu, 
is to determine the distribution parameter from a fit of the lepton spectrum away from the 
endpoint where the spectrum is dominated by 6 ^ c transitions. This result is then used 
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to model the endpoint region, which purely originates from h ^ u transitions, in order to 
extract the value of from the data. Considering however the fact that the decay 

b ^ c involves a quark mass rric, which cannot be neglected, the corresponding distribution 
parameter might be unsuitable to describe simultaneously the transition h ^ u. 

On the other hand, the independent determination of the heavy quark momentum 
distribution from the entire J / ip spectrum of the B ^ J/ip X decay involves the effective 
transition 6 — > s (and also the Cabibbo suppressed decay b ^ d), i.e., a small (current) mass 
for the final state. Therefore the corresponding distribution parameter is still appropriate 
for the extraction of iV^ufe/Kbl, which requires the distribution function associated with the 
transition to a massless final state. 



The CLEO analysis was performed within the ACCMM model using the value 
Pf = 0.3 GeV for both, b c and b u decays. The modification which arises when 
keeping this value for the transition b c but applying Pf = 0.5 GeV for b ^ u reads 





2 




2 




p/=0.5 




p^=0.3 



0.3) 



T{pf = 0.5) 



(40) 



where f (p/) = dEi-^^{pf) denotes the integration over the endpoint domain of the lep- 
tonic spectrum supposing \ Vub\ = 1- This is exactly the same relation which has been used 
in Ref. and we quote the corresponding value for \ Vub/Vcb\'^ in Section |3]^. Nevertheless 
the physical interpretation is somewhat different. The authors of Ref. |3^ took the value 



of \Vcb\ as determined independently from other analyses. This appeared to be necessary, 
because they did not consider the limitations to the universality of the distribution func- 
tion arising from different masses of the final state quark. In contrast to this we argue that 
Eq. (^) can be applied directly in the investigation of the semileptonic spectrum, since 
only the distribution parameter governing the endpoint region has to be changed relative 



to the existing CLEO analysis Within this analysis the application of the ACCMM 
model to 6 — > c transitions has to be regarded as a phenomenological fit, because from a 
theoretical point of view the neglect of Fermi motion in the final state is only appropriate. 



when small final state quark masses are involved |^ . 

The latter remark illustrates the importance of a study which is independent of the 
semileptonic data. Our investigation of the decay B J/tpX allows a suitable determi- 
nation of the parameters of the ACCMM model which can be used for the analysis of the 
endpoint spectrum in the decay B Xlv. The same statement holds true for our analysis 
within the parton model. 
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5 Polarization of the J/^ 



Applying the same method as for the calculation of the unpolarized J/ip spectrum, we 
proceed to determine the momentum spectrum of a longitudinally polarized state. 

In references and the polarization of the J/ip was obtained from a study of 
the free quark decay h — J/ips, i.e., for fixed momenta in the final state. The result, 
r^/r ~ 0.54, was identified with the average polarization in the corresponding decay of the 
B meson. Considering the bound state structure of the B within our inclusive approach 
yields the momentum distribution of the J/ip. Consequently, it allows us to determine 
the polarization in various kinematic regions, where the results can be compared with 
measurements from CLEO fHSU and ARGUS K . 



Since our model is based on local quark-hadron duality, we do not expect to reproduce 
the polarization within the region \k^\ > 1.4 GeV, because it is governed by the two-body 
decay modes B — > J/ipK^*\ As can be concluded from the ARGUS data (see Tab. 1), 
after subtraction of the exclusive mode B ^ J / ipK the high momentum range of the decay 
spectrum is dominated by a single orbital angular momentum and consequently by a single 
CP eigenstate of the mode B —>■ J/ipK* ||^. This feature is not accounted for within our 
inclusive approach since final state interactions are not part of our consideration. 

Nevertheless the study of the low momentum range, which involves nonresonant multi- 
particle states, provides an additional test of our approach, especially in view of future 
measurements with reduced experimental errors. Furthermore we investigate the modifi- 
cation of the average polarization relative to the free quark decay due to the bound state 
structure of the -B, in which the hadronic vertex acts on the polarization of the J /ip through 
the momentum distribution of the underlying h quark. 

Replacing the polarization sum in Eq. (|^) by e^{\ = 0)£:i^(A = O), a calculation analogous to 
the determination of the unpolarized B ^ J/ip X decay spectrum within the parton model, 
introduced in Section |^, yields the momentum spectrum for the longitudinally polarized 
J/ip. In the restframe of the B, 



1 dV 




(41) 



'1 + — U2- ^+ - 7727 V ] ] + {x+^X- 



with \Cf \ defined according to Eq. (18). Note that the difference to the unpolarized spec- 



trum occurs in the first term in the parenthesis of Eq. (0), whose sign changes compared 
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to Eq. pop. The partial branching ratio in the CLEO frame is obtained by performing the 
integration due to the transformation of the spectrum to that for a B meson in flight (see 
Eq. (P^D), as well as the integration over the momentum range which is considered. 

The calculation in the ACCMM model follows the lines of Section ^, where for the 
polarized case Eq. ( ^9]) is replaced by the corresponding expression for the decay width of 
b —>■ {J Qf- In the restframe of the heavy quark one obtains 



<f) 

0,L 



271 



(w 



m 



f 



Ml 



(42) 



The partial branching ratio for the momentum range ki < \k.A < k2 then reads 



(43) 



where 



tb j d\k^\ j dp 



SbJ) 



{p)-\k^-'^\p)\ ^ k^{\k'^\)~\k_{\k'^\ 



d\kL 



92{\ky,\, p) 



{\k-{\h\)l \k^-'^\p)\]K\k^lp) 
in { A;+(|^^|) , kt'^\p) ]h{\k^\,p) 



max 
min 



and 



h{\k^lp) = e k^l\\k^\)-\k'^^'\p)\ 9 kf^\p)-\k.{\k^ 



(bj), 



Ab,f). 



(44) 
(45) 

(46) 



In Eqs. ( ^^ - |46|) the limits of integration are written in a way as to take into account 
simultaneously the Lorentz boost from the h to the B restframe and the subsequent boost 
to the CLEO frame. 

Our results for the average longitudinal polarization of the J/ip., Tl/T = ABl/AB, 
are presented in Tab. 1. There we give the polarization obtained in the parton and the 
ACCMM model for different values of the quark momentum distribution parameters Ep 
and pf, respectively. The J/ip momentum ranges that we considered are chosen according 
to the experimental data from ARGUS and CLEO. From the table one can read that 
in both models the data is reproduced correctly for the range \k^\ < 1.4 GeV and also 
for the average over the complete momentum range. The agreement is good, but the 
experimental errors are still large. On the other hand, the polarization predicted for the 
high momentum region underestimates the data. This is to be expected because, as we 
mentioned, the region > 1.4 GeV is dominated by two-body modes. 
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J/ip momentum 


CLEO II [ 


43] 


ARGUS 


44 




< 0.8 GeV 
0.8 GeV <k^< 1.4 GeV 
1.4 GeV <k^< 2.0 GeV 
all k^ < 2.0 GeV 


0.55 ±0.35 
0.49 ±0.32 
0.78 ±0.17 
0.59 ±0.15 


1.17±0.17 



parton model 


J/ip momentum 


Bp = 0.004 


Ep = 0.006 


Ep = 0.008 


Ep = 0.010 


k^ < 0.8 GeV 


0.416 


0.415 


0.414 


0.413 


0.8 GeV< k^ < 1.4 GeV 


0.520 


0.515 


0.512 


0.508 


1.4 GeV< k^ < 2.0 GeV 


0.557 


0.552 


0.547 


0.543 


all k^ < 2.0 GeV 


0.537 


0.529 


0.522 


0.516 




ACCMM model 


J/ip momentum 


Pf = 0.3 


Pf = 0.4 


Pf = 0.5 


Pf = 0.55 


k^ < 0.8 GeV 


0.515 


0.504 


0.495 


0.491 


0.8 GeV< k^ < 1.4 GeV 


0.548 


0.538 


0.530 


0.527 


1.4 GeV< k^ < 2.0 GeV 


0.556 


0.549 


0.542 


0.539 


all k^ < 2.0 GeV 


0.553 


0.543 


0.534 


0.529 



Table 1: J/ip polarization F^/r in the decayS J/ipX in the parton and the ACCMM 
model compared to data. The parameters are nigp = 0.2 GeV, rris = 0.125 GeV, rrid = 0, 
and Ep, pf[GeV] respectively, as given above. 

Significant differences between the parton and the ACCMM model show up only for 
the low momentum range {\k^\ < 0.8 GeV), which may provide an additional feature for 
distinguishing between the models. The bound state corrections to the average polarization 
of the J/ip (for all < 2.0 GeV) are marginal; a property which is reflected in the weak 
dependence of the polarization on the distribution parameters. As a consequence of the 
cancellation of binding effects in the ratio AB^/AB, the determination of Ep and pf, 
respectively, from the polarization turns out to be impossible. 

We conclude that the inclusive approach to i? —>■ J/ipX decays in the parton and 
the ACCMM model yields the polarization of the J/ip in various kinematic regions. Our 
analysis shows that for the high momentum range, in which flnal state interactions are 
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important, the approach is less rehable in the polarized than in the unpolarized case. 
On the other hand, the inclusive description is well suited in the low momentum range 
governed by nonresonant multi-particle states. A study of the polarization in this range 
may help distinguishing between the models considered, as soon as the experimental error 
is reduced. 

6 Summary and Conclusion 

We incorporate in this article the bound state effects of the B meson in the analysis 
of the direct decay B — > J/ipX. An inclusive approach has been worked out in detail 
within the framework of the parton and the ACCMM model where in each case a one- 
parameter momentum distribution function for the heavy quark is introduced. We fixed 
the distribution parameter by comparing the predicted J/ip momentum spectrum with 
the recent CLEO data, putting emphasis on the adequate reproduction of the sizable low 
momentum spectrum, which contains nonresonant multi-particle final states. Both models 
yield a b quark momentum distribution softer than the one commonly used in studies of 
semileptonic B decays. Within the parton model we obtain Sp — (9(0.008) where moderate 
deviations from the Peterson et al. distribution, taken from production experiments, are 
apparent in the data. In further studies of inclusive B decays within the parton model it 
will be of interest to consider a modified Ansatz for the h quark momentum distribution. 

The ACCMM model can account for the data over the whole range of phase-space when 
we use a Fermi motion parameter pf = 0(0.5 GeV). Applying this result to the cndpoint 
electron energy spectrum of semileptonic B meson decays yields a value 10^ x \Vub/Vcb\'^ — 
1.03, which is in accordance to the result obtained within the exclusive ISGW model. 

The successful reproduction of the experimental data confirms the validity of the fac- 
torization assumption for inclusive B decays where the bound state effects in both models 
imply a large value for the effective color singlet coefficient, \a2\ = C(0.28) (when choosing 
tb = 1.54 ps, \Vcb\ = 0.043). The inclusive approach can also account for the measured av- 
erage polarization of the J/ip which is independent of the normalization constants (|a2|/^) 
and Tb- 

In view of the fact that semileptonic decay spectra away from the endpoint of the 
phase-space involve 6 — > c transitions a study of the decay B ^ J/ip X presently provides 
the only possibility to extract the momentum distribution of the b quark corresponding to 
a decay in which the mass of the quark in the final state is negligible. This extraction is 
direct as no integration over the distribution function has to be performed. Furthermore 
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the momentum spectrum of the J/%1) is more sensitive to the bound state structure of the 
B than the electron energy spectrum from semileptonic decays. The difficulties which arise 
in the study of a momentum spectrum containing resonance structures can be avoided as 
soon as precise measurements are available which allow a determination of the distribution 
function from the endpoint domain of the electron energy spectrum in semileptonic decays 
or the photon spectrum in 6 — > s 7 decays. 
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Figure Captions 

Fig. la Dominant contribution to the decay B J/ijj X with the light final state quark 
(/ = s, d) in the subprocess having positive energy. 

Fig. lb Contribution to the decay S — > J/il^ X with the light final state quark qf (/ = s, d) 
in the subprocess having negative energy. 

Fig. Ic b ^ J/ip + s{d) matrix element. The square represents a four-fermion vertex. 

Fig. 2a Momentum (respectively mass) distribution function f£{x) in the IPM for various 
values of Sp. 

Fig. 2b Mass distribution function w{x) in the ACCMM model for various values of p/. 
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Fig. 3 x±{f — s) as functions of the J/tp momentum for = 0.125 GeV (current mass) 
and for 777.5 — 0.55 GeV (constituent mass). 

Fig. 4 Theoretical momentum spectrum in the IPM for direct inclusive J/ ip production from 
B decays at the T(4s) resonance. The vahies for the parameters are = 0.125 GeV, 
md = (current masses), Mb = 5.279 GeV, = 0.384 GeV, jaal = 0.275, tb = 1.54 
ps, |V"c6| = 0.043, \Vcs\ = 0.97, \Vcd\ = 0.22 and various values of Bp as shown. 

Fig. 5 Same as in Fig. 4, now theoretical momentum spectrum for various values of Ep and 
1 02! as shown, compared with the CLEG data. 

Fig. 6 Same as in Fig. 5, now for Sp = 0.008, \a2\ = 0.285 and current, respectively con- 
stituent masses for the final state quarks. 

Fig. 7 Theoretical momentum spectrum in the IPM for direct inclusive J/ip production 
compared with the corrected GLEO data where the expected contribution of the 
exclusive two-body final states {J/ipK'^*^) has been subtracted. The values for the 
parameters are 777^ = 0.55 GeV, = 0.33 GeV, |a2| = 0.275 and Sp as shown. 

Fig. 8 Theoretical momentum spectrum in the AGGMM model for direct inclusive J/ip 
production from B decays at the T(4s) resonance compared with the GLEO data. 
The values for the parameters are 777^ = 0.125 GeV, 777^ = (current masses), nisp — 
0.2 GeV, Mb = 5.279 GeV, tb = 1.54 ps, \Vcb\ = 0.043, |V;,| = 0.97, \Vcd\ = 0.22 and 
various sets of values for pf and \a2\ as shown. 

Fig. 9 Same as in Fig. 8 , now for pf = 0.55 GeV, \a2\ = 0.28, rusp — 0.2 GeV and current, 
respectively constituent masses for the final state quarks. 

Fig. 10 Theoretical momentum spectrum for direct inclusive J/ip production in the AGGMM 
model for pf — 0.55 GeV, \a2\ — 0.28 and various values of shown compared 

with the GLEO data. 

Fig. 11 Theoretical momentum spectrum as in Fig. 10, now for TTigp = 0.2 GeV, \a2\ = 0.28 
and various values of p/ as shown. 
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fig. la 



fig. lb 



fig. Ic 
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